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A search for direct production of scalar bottom quarks (6) is performed with 310 pb _1 of data 
collected by the D0 experiment in pp collisions at = 1.96 TeV at the Fermilab Tevatron Collider. 
The topology analyzed consists of two b jets and an imbalance in transverse momentum due to 
undetected neutralinos (x?), with x? assumed to be the lightest supersymmetric particle. We find 
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the data consistent with standard model expectations, and set a 95% C.L. exclusion domain in the 
(mj,trijo) mass plane, improving significantly upon the results from Run I of the Tevatron. 



PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 

Supersymmetric (SUSY) models provide an exten- 
sion of the standard model (SM) with mechanisms viable 
for the unification of interactions and a solution to the hi- 
erarchy problem. Particularly attractive are models that 
conserve imparity, in which SUSY particles are produced 
in pairs and the lightest supersymmetric particle (LSP) 
is stable. In SUSY, a scalar field is associated to each of 
the left and right handed chirality states of a given SM 
quark or lepton. Two mass eigenstates result from the 
mixing of these scalar fields. The spin-1/2 partners of the 
neutral gauge and Higgs bosons are called neutralinos. 

In supergravity inspired models Q, the lighest neu- 
tralino Xi arises as the natural LSP, and, being neutral 
and weakly interacting, could be responsible for the dark 
matter in the universe. For large values of tan /3 (the ratio 
of the vacuum expectation values of the two Higgs fields) 
the mixing term among the scalar fields associated with 
the bottom quark is large. Therefore, a large splitting 
is expected among the mass eigenstates, that could re- 
sult in a low mass value for one of them, hereafter called 
scalar bottom quark or sbottom (b). The SUSY particle 
mass hierarchy can even be such that the decay b —>■ b\\ 
is the only one kinematically allowed Q, an assumption 
that is made in the following. 

In this Letter, a search is reported for b pair production 
with 310 pb" 1 of data collected during Run II of the Fer- 
milab Tevatron. At leading order, the b pair production 
cross section in pp collisions depends only on the sbottom 
mass. For a center-of-mass energy y/s = 1.96 TeV, the 
next-to-leading order (NLO) cross section, calculated with 
PROSPINO-2 jj| ranges from 15 to 0.084 pb for sbottom 
masses between 100 and 230 GeV, with very little depen- 
dence on the masses of the other SUSY particles. The 
final state of this process corresponds to two b jets and 
missing transverse energy due to the undetected 

neutralinos. The maximum sbottom mass (m^) excluded 
by previous results is 148 GeV ||. 

A full description of the D0 detector is available in 
Ref . . The central tracking system consists of a silicon 
microstrip tracker and a central fiber tracker, both lo- 
cated within a 1.9 T superconducting solenoid. A liquid- 
argon and uranium calorimeter covers pseudorapidities 
up to \r)\ ~ 4.2, where rj = — In [tan (0/2)] and 9 is the 
polar angle relative to the proton beam. The calorimeter 
has three sections, housed in separate cryostats: the cen- 
tral one covers |?7| < 1.1, and the two end sections extend 
the coverage to larger \rj\. The calorimeter is segmented 
in depth, with four electromagnetic layers followed by up 
to five hadronic layers. It is also segmented into pro- 
jective towers of 0.1 x 0.1 size in rj — <f> space, where 4> 



is the azimuth in radians. An outer muon system, cov- 
ering |?7| < 2, consists of a layer of tracking detectors 
and scintillation trigger counters positioned in front of 
1.8 T toroids, followed by two similar layers after the 
toroids. Jet reconstruction is based on the Run II cone 
algorithm Q with a cone size of 0.5, that uses ener- 
gies deposited in calorimeter towers. Jet energies are 
calibrated using transverse momentum balance in pho- 
ton+jet events. The missing transverse energy in an 
event is based on all calorimeter cells, and is corrected for 
the jet energy calibration and for reconstructed muons. 

The D0 trigger has three levels: LI, L2, and L3. The 
data were collected with triggers specifically designed for 
^ T +jets topologies. We define ftp = l^jets^l tne 
vector sum of the jet transverse momenta. The trig- 
ger conditions at LI require at least three calorimeter 
towers with Et > 5 GeV, where a trigger tower spans 
x Arj = 0.2 x 0.2. We then require ftj. > 20 (30) GeV 
at L2 (L3). Approximately 14 million events were col- 
lected with the $ T +jets triggers. 

The signal is simulated in the framework of a generic 
minimal supersymmetric standard model, in which we 
vary the masses of the b and Xii au other parameters be- 
ing fixed. The masses of the other SUSY particles are set 
such that the only sbottom decay mode is into b\i- The 
SUSY and SM processes are processed using Monte Carlo 
(MC) generators PYTHIA 6.202 8] for the signal, ALP- 
GEN 1.3.3 interfaced with pythia for the SM. All the 
events are passed through a full GEANT-3 10] simulation 
of the geometry and response of the D0 detector with 
an average of 0.8 minimum-bias events overlayed on each 
generated event. The CTEQ5L parton density functions 
(PDF) [13 are used in the simulation. 

Instrumental background from mismeasurement of jet 
energies in multijet events is estimated from data, and 
is referred to as "QCD" background in the following. The 
main SM backgrounds relevant to our analysis are from 
vector boson production in association with jets, and top 
quark production. To estimate the backgrounds from 
W/Z+jets processes, we use the NLO cross sections com- 
puted with MCFM 0. The theoretical NNL0 ti produc- 
tion cross section is taken from Ref. . 

The events are examined to ensure that the recon- 
structed vertex corresponds to the actual position of the 
primary vertex (PV). We select events that are well con- 
tained in the detector by restricting the PV within 60 cm 
along the beam direction with respect to the detector 
center. We define a charged-particle fraction (CPF) as 
the ratio of the charged-particle transverse energy, com- 
puted from the sum of scalar px values of charged parti- 
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cles (reconstructed in the tracking system) that emanate 
from the PV and are associated with a jet, divided by 
the jet transverse energy measured in the calorimeter. 
The two leading jets, i.e those with the largest trans- 
verse energies, are required to have CPF>0.05. This cri- 
terion rejects events with fake jets or where a wrong PV 
is selected. The overall inefficiency associated with this 
procedure is measured using events collected at random 
beam crossings, and events with two jets emitted back-to- 
back in azimuth. The jets must also have energy fraction 
in the electromagnetic layers of the calorimeter <0.95 
and pt> 30, 15 GeV for the first and second leading jets. 
This set of initial cuts requires in addition A<f> < 165°, 
where A<p is defined as the difference in azimuth between 
the two leading jets. 

Table H| defines our selection criteria, and shows the ef- 
fect of applying them sequentially in the analysis of data, 
and their impact on signal efficiency, for the choice of 
(mg,m^o) = (140,80) GeV. Criteria C1-C4 are effective 
against QCD, C2 and C4-C8 against vector-bosons+jets, 
and C9 suppresses tt background. For 6 masses of 
~ 100 GeV, the mean $ T and jet px are close to what 
is expected from SM backgrounds, but are substantially 
larger for higher 6 masses. The selections are tuned on 
MC so as to maintain good sensitivity to signal for small 6 
masses, using minimal values for threshold requirements, 
for instance # T >60 GeV (CI) and p T > 40, 20 GeV (C2) 
for the first and second leading jets. Later we show that, 
depending on the masses (mg,m^o), higher thresholds on 
frp and jet px can be applied to increase the sensitivity 
to signal. 

The first and second leading jets are required to be in 
the central region of the calorimeter, |r? det | <1.1 and 2.0 
respectively (C3), where ?7 dct is the jet pseudorapid- 
ity calculated with a jet origin at the detector cen- 
ter. Because of the central production of b events, 
these selections do not affect signal efficiency, but re- 
duce background. We also define A0 m ,- n (^ r ,jets) and 
A<j) max {$ T ,\ets) , the minimum and maximum of the dif- 
ferences in azimuth between the direction of $ T and the 
direction of any jet. Requiring A(j) min >35° rejects QCD 
events (C4), and A(j) min <120° and A(j) max <175° sup- 
press SM background (C4, C5). 

Since we do not expect isolated electrons, muons or 
tau leptons in signal events, vetoes are imposed on events 
with an isolated electron (C6), muon (C7), or a charged 
track (C8) with px>5 GeV. Electrons and muons are de- 
fined isolated based on a criterion for energy deposition 
in a cone around the lepton direction in the calorimeter. 
A charged track is considered isolated if no other charged 
track with p-r>1.5 GeV is found in a hollow cone with 
inner and outer radii 0.05 and 0.2, formed around the 
direction of the track. The last requirement (C9) stipu- 
lates that either two or three jets are allowed. 

Table [n] gives the numbers of events expected for SM 
backgrounds and signal, and the number of events ob- 



TABLE I: Sequence of criteria applied for the selection of 
events with their corresponding impact on data and on signal 
efficiency (Eff.) for (m^m^o) = (140,80). 



Selection criterion applied Events left Eff. (%) 



CI 


$ T > 60 GeV 


16,279 


18 


C2 


p Tl > 40 GeV, pr 2 > 20 GeV 


14,095 


16 


C3 


IvfctW < 1.1, fog&l <2.0 


9,653 


14 


C4 


35° < A(t> m in0 T , jets) < 120° 


3,149 


10 


C5 


Ac/We (# T , jets) < 175° 


2,783 


9 


C6 


isolated electron veto 


2,059 


9 


C7 


isolated muon veto 


1,809 


9 


C8 


isolated track veto 


756 


7 


C9 


2 or 3 jets 


671 


6 



served in data after the above selections. Since an impor- 
tant fraction of the background corresponds to processes 
with light-flavor jets in the final state, we take advan- 
tage of the presence of b jets in the signal to significantly 
increase the sensitivity of the search by using a lifetime- 
based heavy-flavor tagging algorithm (6-tagging). Based 
on the impact parameters of the tracks in the jet, the al- 
gorithm [14j computes a probability for a jet to be light- 
flavored. 

We select the 6-tagging probability such that 0.1% of 
the light-flavored jets are tagged for jets having p? of 
50 GeV as yielding the best expected signal sensitivity. 
The corresponding typical tagging efficiencies for c- and 
6-quark jets are 5% and 30%, respectively. Because the 
current detector simulation does not reproduce the track- 
ing precisely enough, the b-tagging algorithm is not ap- 
plied to simulated jets directly. Instead, jets are weighted 
by their probability to be 6-tagged, according to their fla- 
vor, using parameterizations derived from data. In what 
follows, we require at least one 6-tagged jet in the event. 
Requiring more than one 6-tagged jet would lower slightly 
the sensitivity of the analysis. 

In order to estimate the background from QCD, we com- 
pare our selected data sample, without imposing the cri- 
terion on Ifjrp (CI), to the simulation of background from 
SM. Figure shows that data are well reproduced by the 
SM background at high $J T . We therefore attribute the 
exponential rise at low l£ T to QCD multijet instrumen- 
tal background. A fit by an exponential to the data for 
lp T < 60 GeV, after subtraction of the contributions from 
the SM, is shown in the insert in Fig. ^ When the fit is 
extrapolated to $J T > 60 GeV, it provides an estimate of 
109±9 QCD events. After 6-tagging, this procedure esti- 
mates the presence of only 4±2 events. Given the larger 
]fi T threshold we use for higher sbottom masses, we ex- 
pect that, after the 6-tagging, less than two QCD events 
will survive the final event selection. The QCD contri- 
bution is therefore neglected in the rest of this analysis. 
Table ITT1 shows the results after all selections, including 
6-tagging, for SM backgrounds, data and signal. 

As already mentioned, the mean $ T and jet pr become 
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TABLE II: Numbers of events expected from SM and QCD 
backgrounds, of data events observed, and of signal events 
expected, after all selection criteria, both before (N exp ) and 
after 6-tagging. All uncertainties are statistical only. Back- 
grounds from 6, c, and light jets (j) are shown separately. 



SM process 




exp 


with 6-tagging 


W(ev + fiv) + jj 




155 ± 13 


1.9 ± 0.2 


W(ev ± fiv) ± cj 




2.2 ± 0.6 


0.2 ± 0.1 


W(ev + fj,u) ±66 




1.1 ± 0.1 


0.6 ± 0.1 


W(tv)+ > 1 jet 




101 ± 14 


4.1 ± 0.6 


W{tv) ± bb 




2.2 ± 0.3 


1.0 ± 0.1 


Z(uP)+jj 




257 ± 12 


3.9 ± 0.2 


Z(yv) ± cc 




8.0 ± 0.7 


0.9 ± 0.1 


Z{vv) + bb 




7.8 ± 0.3 


4.0 ± 0.2 


WW, WZ, ZZ 




14.2 ± 0.7 


0.9 ± 0.2 


top production 




7.9 ± 0.2 


3.8 ± 0.2 


Total SM 




556 ± 23 


21.5 ± 0.8 


QCD background 




109 ± 9 


4 ± 2 


Data 




671 


22 


(m s ,m s o) = (140, 8( 


)) GeV 


43 ± 2 


23.1 ± 0.9 
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FIG. 1: Distribution in lp T after applying all criteria, except 
$ T > 60 GeV (CI). The dark shaded area corresponds to the 
SM simulation. A fit by an exponential to $ T < 60 GeV, after 
subtraction of the contributions from the SM, shown in the 
figure insert, is used to estimate the instrumental background. 



substantially larger for higher sbottom masses than the 
values expected from SM backgrounds. This provides a 
handle for improving the sensitivity to the signal for large 
mg. Table ILTTI shows results for two higher sbottom-mass 
points, the chosen $ T and px thresholds, together with 
the resulting number of events found after all selections, 
including 6-tagging, for data, SM background and signal. 
For the highest sbottom masses probed, we note a deficit 
in the number of events observed compared to the SM 
background expectation. The probability of such a deficit 
is 4%. 

The following systematic uncertainties are taken into 
account in deriving the final results. The integrated lu- 



TABLE III: Optimized values for the criteria CI and C2, 
numbers of data events observed, numbers of events expected 
from SM and signal for two (m^m^o) masses after 6-tagging 
(statistical uncertainties only). 



(ra^m^o) in GeV 


(180,90) 


(215,0) 


CI: # T [GeV] 


60 


80 


C2: jet 1 p T [GeV] 


70 


100 


C2: jet 2 p T [GeV] 


40 


50 


data 


7 





SM 


8.9 ±0.3 


3.2 ±0.2 


signal 


9.4 ±0.3 


4.6 ±0.1 



minosity contributes an uncertainty of 6.5%. The uncer- 
tainty from jet energy calibration is typically of the order 
of 7%. The total uncertainty from jet energy resolution, 
jet track confirmation, misvertexing and jet reconstruc- 
tion is 5%. The systematic uncertainties from NL0 cross 
sections in the SM backgrounds are estimated to be 15%. 
The effect of the choice of PDF on signal efficiencies is 
evaluated using the CTEQ6 . 1M PDF error set 0] resulting 
in a 8% uncertainty. The uncertainty from MC statistics 
can reach 10% for the SM and 5% for signal. The total 
uncertainty from isolated electron, muon, and track ve- 
toes is 9%. The uncertainty from heavy- flavor tagging is 
12% for SM and 8% for signal. Finally, the uncertainty 
from the trigger efficiency is 5%. 

Since we do not observe any excess in the data rel- 
ative to the expectations from SM backgrounds, we set 
limits on the production of sbottom quarks. Observed 
and expected 95% confidence level (C.L.) cross section 
upper limits are obtained using the modified frequentist 
approach |l6| . with correlations included between sys- 
tematic uncertainties. The NL0 6 pair production cross 
section is subject to theoretical uncertainties arising from 
the PDF and from the renormalization and factorization 
scale choices. For a b mass of 200 GeV, a 16% PDF un- 
certainty is evaluated using the CTEQ6 . 1M PDF error set, 
and a 12% uncertainty is found by varying the scale by 
a factor of two up or down. For a given neutralino mass, 
a sbottom mass limit is obtained where the cross sec- 
tion upper limit intersects the production cross section 
reduced by these uncertainties combined in quadrature. 
The results are summarized in the 95% C.L. exclusion 
contours displayed in Fig. [3 At higher sbottom masses, 
no events are observed where about three are expected, 
leading to an observed limit more constraining than ex- 
pected. 

In summary, this analysis represents the first Tevatron 
Run II search for pair production of scalar bottom 
quarks. The exclusion contour we obtain is substantially 
more restrictive than the ones published with Run 
I Tevatron data. With the current analysis using 
310 pb _1 , the maximum excluded is 222 GeV, an 
improvement of more than 70 GeV with respect to 
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FIG. 2: Excluded regions at the 95% C.L. in the sbottom and 
neutralino mass plane. The new region excluded by this anal- 
ysis is shown in dark shading. The dashed line corresponds to 
the expected limit. Regions excluded by previous experiments 
are also displayed in the figure |f|- 



previous results, and the most restrictive limit on the 
sbottom mass to date. 
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